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Abstract: A study was conducted to investigate the land-use pattern change over a period of 18 years (1988-2006) by using remote sens¬ 
ing and Geographical Information System (GIS) technologies, in the West Bhanugach Reserved Forest, a hill forest, in Sylhet Forest Divi¬ 
sion of Bangladesh. The images were processed using ERDAS Imagine software. Both supervised and unsupervised approaches were ap¬ 
plied and ground control points were collected using a GPS. Maps were prepared using GIS software. Results showed that vegetation 
cover drastically decreased from the year 1988 to 1996 (1 826 ha to 1 714.85 ha), but increased gradually from the year 1996 to 2006 
(1 714.85 ha to 1 847.83 ha) due to the initiation of co-management practice involving local communities. Change in bare land was in¬ 
versely proportionate to the amount of vegetation cover changes unless any other land-uses were converted into bare land. The area of 
water bodies increased from the year 1988 to 1996 (307.67 ha to 379.53 ha), but decreased from the year 1996 to 1997, then remained in- 
variabile from the year 1997 to 2006. Some recommendations were also made for applying the RS and GIS techniques to study the 
land-use pattern change in the Bhanugach Reserved Forest and to create a GIS data base for the study area. 
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Introduction 

The landscape is dynamic in relation to spatial, structural and 
functional patterns. Land-use is the purpose for which human 
exploit the land cover (Fresco 1994). It is a product of interac¬ 
tions between a society's cultural background, state, and its 
physical needs on the one hand, and the natural potential of land 
on the other (Ram et al. 1993). Resources, ecosystem, biophysi¬ 
cal environment, and land-use on the surface of the earth un¬ 
dergo changes over time. It is obviously constrained by envi¬ 
ronmental factors (e.g. soil characteristics, climate, topography, 
and vegetation) (McCracker et al. 1998) and human forces (Fer¬ 
nandez et al. 1992) (e.g. specific management practices, social, 
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political and economical forces) (Medley et al. 1995). 

Land-use change is the modification in the purpose of the land, 
which is not necessarily only the change in land cover but also 
changes in intensity and management (Verburg 2000). Both hu¬ 
man-induced and natural land cover/use changes are critical is¬ 
sues due to their great influence on global warming, loss of bio¬ 
diversity, and impact on human life (Sala et al. 2000; Geoghegan 
et al. 2001; Lopez et al. 2001). Change detection is a process of 
identifying and analyzing the differences of an object or a phe¬ 
nomenon through monitoring at different times (Singh 1989; 
Mouat et al. 1993). A wide range of applications can be bene¬ 
fited from the study of change process over a specified area at 
different times. Recent literature deals with the application of 
change analysis to different problem domains (e.g. land 
use/cover change dynamics, global change analysis, monitoring 
of pressure on the environment, monitoring of agricultural pro¬ 
duction, assessing damages due to forest fires and deforestation, 
and monitoring damages due to natural calamities like floods, 
earthquakes, and volcanic eruption etc.) (Anwar 2002). Informa¬ 
tion about land-use change is necessary to update land cover 
maps and for effective management and planning of the re¬ 
sources for sustainable development (Alphan 2003; Muttitanon 
et al. 2005). Furthermore, timely and reliable data on land-use 
may facilitate the formation of integrated resource management 
policies. Even though detailed land-use information are well 
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equipped and updated in most of the developed countries, lack of 
geospatial database for land-use persists in developing countries 
and may hinder proper development planning. 

Bangladesh is a developing country experiencing high popula¬ 
tion growth where unemployment has been exacerbating the 
overall situation with a great pressure on natural resources. For¬ 
est is an important one among other natural resources but due to 
extreme poverty the forests there are under the threat of defores¬ 
tation (BBS 2004). In hilly region of the country, this scenario is 
ruder (Rasheed 2003). Because, the hill peoples have very lim¬ 
ited alternative sources of income but depend on forests. Due to 
their inconsistent use (shifting cultivation, intentional fire, over 
exploitation, illegal felling etc.) of the forests they are making 
the unsustainable use of these natural resources (Iftekhar 2006). 
According to Allan (1986), hilly areas are particularly sensitive 
to drivers of land-use changes. Land abandonment and forest 
expansion is common in hilly areas worldwide. But in case of 
Bangladesh, due to high population pressure and other 
socio-cultural reasons, the hill peoples will not abandon their 
places creating no scope for forest expansion; rather, with the 
hunger-stricken situation, they along with the migrated mainland 
people do harm to forests regularly. In this circumstance, the idea 
of sustainable development seems to be abstract with the avail¬ 
ability of no or very few updated information on land-use pattern 
and changes in the hilly region of the country. For sustainable 
utilization of the forested ecosystems, it is essential to know the 
natural characteristics, extent and location, its quality, productiv¬ 
ity, suitability and limitations of various land-uses. In order to 
improve the economic condition of the area without further dete¬ 
riorating the bio-environment, every bit of the available land has 
to be used in the most rational way. This requires present and 
past landuse/land cover data of the area (Chaurasia et al. 1996). 

Besides, the economic and ecological importance of forests is 
well understood (Franklin et al. 2000). Thus the importance of 
mapping, quantifying, and monitoring changes in the physical 
characteristics of forest cover has been widely recognized as a 
key element in the study of global change (Nemani et al. 1996). 
Alongside, the location and rates of forest structural change and 
the degree to which landscapes respond to these disturbances 
require extensive investigation (Borak et al. 2000; Rogan et al. 
2002 ). 

In recent years, a growing number of studies have combined 
principles from landscape ecology with spatial-econometric 
methods to account for how human decision-making, ecosystem 
function, and their interaction affect landscape changes across 
different spatial scales (e.g. Kline et al. 2001; Irwin et al. 2002; 
Anwar 2002; Hite et al. 2003; Fang et al. 2005; Srinivisan 2005; 
Vance et al. 2005), many of which have used data from satellite 
imagery to study the land-use-change pattern, and information 
obtained from these studies can illuminate the still poorly under¬ 
stood role of natural amenities as drivers of development (Alig et 
al. 2004). 

Application of remotely sensed data made possible to study 
the changes in land cover in less time, at low cost and with better 
accuracy (Kachhwaha 1985) in association with Geographical 
Information System (GIS) that provide suitable platform for data 
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analysis, update and retrieval (Star et al. 1997; McCracker et al. 
1998; Chilar 2000). Hence, space-bome remotely sensed data 
may be particularly useful in developing countries where recent 
and reliable spatial information is lacking (Dong et al. 1997). In 
Bangladesh a few studies (Islam et al. 1988; Quadir et al. 1998; 
Hossain 2005; Islam et al. 2006) on land-use change have been 
conducted hitherto leaving a greater scope for working on this 
sector for ensuring sustainable development. Considering the 
present demand and urgency, the present study was conducted in 
the West Bhanugach Reserved Forest in Sylhet Forest Division 
of Bangladesh to investigate the land-use pattern change over a 
period of 18 years (February, 1988- February, 2006) using re¬ 
mote sensing and GIS technologies, which will be helpful in 
creating dynamic information system for land-use planning and 
also in developing national development policy. 

Materials and methods 

Study Site 

A deliberate sampling method was applied to select the West 
Bhanugach Reserved Forest (Fig. 1) as the study area. The site is 
located in the civil administrative units of Kamalganj Thana, 
Maulvibazar District under the forest administration unit of 
Moulavibazar range under Sylhet Forest Division, Bangldesh. 
The GPS location of the site is between 24°18 / -24°24' N and 
91 ° 46 '- 91 ° 50 ' E. The reserved forest is divided into three beats; 
viz. Lawachara, Chautali and Kalachara with a total area of 
2 610.31 ha (Anon 1970), of which Lawachara beat was declared 
as a national park in 1996 (Gazette Notification (PBM 
(S-3)7/96/367 on 07 July 1996) with an area of 1 250 ha (NA- 
COM 2003). The other beats, Kalachara and Chautali, are of 
878.2 ha and 482.11 ha, respectively. 

This forest falls under the bio-ecological zone: 9b-Sylhet Hills 
(Nishat et al. 2002). The soils are brown, sandy clay loam to clay 
loam of Pliocene origin (Hossain et al. 1989). The area is undu¬ 
lating with slopes and hillocks ranging from 10-50 m and are 
scattered in the forest. Numerous streams flow through the forest. 
Numerous trails and tracks are found within the forest, created 
by the local people for collecting wood from the forest (Feeroz et 
al. 2000). 

The climate of the forest is in general warm and humid but the 
weather is cool and pleasant during winter. The temperature 
varies on an average from minimum 26.8°C in February to 
36.1°C in June. The humidity is high in the forest throughout the 
year, with monthly average humidity varying from 74% in 
March to 89% in July. The area covered under the reserved forest 
is one of the wettest in the country, thus, the rainfall is quite high, 
with an annual average of 4 000 mm approximately, with maxi¬ 
mum rainfall falling during June to September from South-West 
monsoon. Pre-monsoon Nor’westerly and cyclonic storms are 
accompanied by high-speed winds and rains, which do consid¬ 
erable damage to property and trees (Nishorgo 2007). 

The forest can be classified as semi-evergreen and mixed de¬ 
ciduous forest (Champion et al. 1965; Das 1989; Feeroz et al. 
2000; Ahsan 2000) where tall trees are deciduous and the un- 
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der-storey evergreen (Ahsan 2000). The forest originally sup¬ 
ported an indigenous vegetation cover of mixed tropical ever¬ 
green forest (Alam 1998). However, almost all of the original 
forest cover has been removed or substantially altered and thus 
has turned into a secondary forest. The old plantations (around 
1920’s) are of primarily native species (but many of them are 
exotic to the forest), which have developed a multistoried struc¬ 
ture, including re-growth of creepers and naturally occurring 
trees and undergrowth species. In the oldest of these areas, the 
vegetation cover has demonstrated the structure of natural forest, 
and evolution towards a natural structure (FSP 2000; Chemonics 
2000). The canopy height varies from 10 m to 30 m. The top 
canopy comprises Tectona sp., Ariocarpus chaplasha, Tetrame- 
les sp., Hopea odorata., Toona ciliata, Pygenum sp., etc. The 
second canopy comprises Quercus spp., Syzygium sp., Gmelina 
sp., Dillenia sp., Grewia sp.. Ficus sp., etc. The underneath in¬ 
cludes Bambusa spp., Alsophila sp., Geodorum sp., Eupatorium 
odoratum etc. and several ferns and epiphytes (Feeroz et al. 
2000 ). 

West Vanugach Reserve Forest 


9P4ff • 91'48" 91*50' 



Fig. 1 Showing the study area 


Methods 
Materials used 

Satellite imagery: a) Landsat TM FCC February, 1988 (Band 5, 
3, 4); b) SPOT FCC February, 1996 (XS1, XS2, XS3); c) Land- 
sat TM FCC January, 1997 (Band 2, 3, 4); d) Landsat TM FCC 


February, 2006 (Band 3, 4, 5); 

The imageries have moderate color contrast and are cloud-free. 
Use of Landsat TM and SPOT imageries for vegetation cover 
mapping and change detection is also supported by a number of 
national and international studies: Islam et al. (1983), Naithani 
(1990), Rosenholm (1993), Quadir et al. (1998), Islam et al. 
(2006). 

Forest Cover Maps (scale 1:15 840). The Forest Cover Maps 
are based on the interpretation of aerial photographs from 1981 
and present different forest types. The maps have been prepared 
under the project "Assistance to the Forestry Sector UNDP/FAO 
Project”, BGD/78/017; 

Survey of Bangladesh Topographic Maps (scale 1 : 250,000 
and 1 : 50,000); 

Software used: a) PC Based IMAGINE Ortho BASE software 
(8.3.1 & 8.4) available at SPARRSO laboratory; b) ARC/INFO 
4.0; c) ARCVIEW 3.2 a.; 

GPS: GARMIN GPS 12XL. 

Geometric correction 

To make the images workable, the data needed to be transformed 
to a uniform ground co-ordinate system of a chosen map projec¬ 
tion. With the help of topo maps of 1:50 000, a good number of 
Destination Ground Control Points (GCPs) were selected mainly 
from distinctive features like rail-road crossing, road-river cross¬ 
ing, etc., and the corresponding points were selected in the 
source images. Then the topo-points were digitised and projected 
from geographic to BTM (Bangladesh Transverse Macerator). 
Because of having high resolution (20 m x 20 m) of SPOT it was 
used as reference image for geometric correction of the other 
images. At least 10 GCPs (Ground Control Point) were taken in 
Polynomial 1, and 8 GCPS for up to Polynomial 5. 

Interpretation key 

Interpretation key (Table 1) was developed by using the existing 
maps of the forest cover compared with the tonal characteristics 
of the hardcopy prints/transparencies of the false color composite 
(FCC) imageries. 


Table 1 . Interpretation key for interpretation of Landsat TM and 
SOPT imageries 



Feature 

ID 

Feature 

Signature in the image 


1 

Vegetation 

Dark to light red with 

smooth texture 

TM interpretation 

key 

11 

Bare land 

Bright white with smooth 

texture 


2 

Water 

bodies 

Dark to light blue with 

smooth texture 


1 

Vegetation 

Dark red to bright red with 

smooth texture 

SPOT interpreta- 


11 

Bare land 

Whitish 

tion key 


2 

Water 

bodies 

Dark light blue with smooth 

texture 
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Interpretation of the satellite imageries 

Interpretation of imageries was done using Erdas Imagine with 
the help of available maps and field visits. The images (both 
SPOT and Landsat TM) were brought into the same spatial reso¬ 
lution of 30 m. 


Data analysis 

Data analysis was done using Microsoft Excel 2003 and SPSS 

10 . 

Results and discussions 


Creation of the GIS data base 

The boundaries were digitized from topographic map sheets 
(1:50 000 scale) covering the study area. The imagery interpreta¬ 
tions were digitized, and all the coverages were transferred into 
the real world coordinate system (Lambert conic conformal sys¬ 
tem, Everest Projection). The polygon topologies were built for 
all the above coverages. 

Area measurement and change detection 

From the IMAGINE attribute table, areas from all the maps were 
calculated. The area obtained from the management plan from 
the year 1963-1964 to 1982-1983 of the study site was com¬ 
pared with the area obtained from the Landsat TM imagery of 
the year 1988. It was more or less (deviation not more than 5%) 
similar between them and it was considered as base area for 
change detection. 


The maps were prepared (Fig. 2) from the geometrically cor¬ 
rected images and presented in the Lambert Conic Conformal 
(LCC) (Everest Projection) grid system. While managing the 
polygons in different maps obtained from images of different 
years we made a neighboring. We have merged any polygon less 
than 1 000 m 2 with the nearest one. After measuring the area, the 
areas obtained from the images are matched with the area avail¬ 
able in literatures. The amount of error was 0.4% as measured. 
From the Landsat TM image of the 1988 of the study site, we 
calculated the area as 2 599.08 ha, which is a little bit deviated 
from the area obtained from literature (Anon 1970) as 2 610.31 
ha. The areas of different maps were computed from the vector 
attribute table of the Imagine. The results obtained from the area 
measurement are given in the Table 2. 



Fig. 2 The land-use patterns of Bhanugach Reserved Forest in the year 1988 (a), 1996(b), 1997(c) and 2006 (d), Bangladesh 


Table 2. Showing the areas measured from the maps of different 
land-uses in West Bhanugach Reserved Forest, Bangladesh. 


Year 

Vegetation 

cover (ha) 

Grid code (1) 

Bare land (ha) 

Grid code (11) 

Water bodies 

(ha) 

Grid code (2) 

Un-classed 

(ha) 

Grid code (0) 

Total 

(ha) 

1988 

1826.14 

465.27 

307.67 

NC* 

2599.08 

1996 

1714.85 

504.37 

379.53 

0.33 

2599.08 

1997 

1799.37 

527.78 

271.55 

0.38 

2599.08 

2006 

1847.83 

458.57 

292.37 

0.31 

2599.08 


*NC: Not Calculated. 

As shown in Fig. 3, it was observed that vegetation cover was 
drastically reduced from the year 1988 to 1996. Conversely, it 
had been increased gradually from 1996 to 2006. This is the most 
interesting and hope-inducing finding of the study in that in La- 
wachara National Park, since 1996 co-management practice ini¬ 
tiated by Nishorgo Support Project of BFD has been practiced 
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(Nishorgo 2007). A number of initiatives have been taken to 
motivate local people and encourage them in AIG (Alternative 
Income Generation) activities. As a result, the extent of vegeta¬ 
tion cover has been increased although in a slow rate. 


1900 
a ^1850 
•2 J 1800 
§ ¥ 1750 
° 1700 
1650 
1600 

1988 1996 year 1997 2006 

Fig. 3 The trends in vegetation cover change of West Bhanugach 
Reserved Forest of Bangladesh from the year 1988 to 2006 

Regarding the trends of change in bare land, it was observed 
that the change is inversely proportionate to the amount of vege- 
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tation cover changes unless any other land-uses (e.g. water bod¬ 
ies) are converted into bare land. As the extent of vegetation 
cover increases, the extent of bare land decreases (comparing 
between Fig. 3 and Fig. 4). In the year 1988, the extent of bare 
land was 465.27 ha and it became 504.37 ha in 1996 (Fig. 4). 
Further the extent increased up to 527.78 ha in 1997 due to the 
decrease in water bodies. But, from the year 1997, the extent of 
bare land has drastically been decreased though the extent of 
water bodies remained more or less similar. 



Year 

Fig. 4 The trends of change in bare land of West Bhanugach Re¬ 
served Forest of Bangladesh from the year 1988 to 2006 

The area of water bodies increased from the year 1988 (307.67 
ha) to 1996 (379.53 ha) but decreased from 1996 to 1997. From 
the year 1997 to 2006, the extent of water bodies remained more 
or less the same (Fig. 5). Therefore, it can be inferred that vege¬ 
tation cover of the study site has been increasing with time, 
which reflects the efficiency of co-management approach in 
reserved forest and protected area management in Bangladesh. 



Year 

Fig. 5 The trends of change in water bodies of West Bhanugach Re¬ 
served Forest of Bangladesh from the year 1988 to 2006 

Conclusions 

It is greatly feasible to use RS and GIS techniques to study the 
changes in land-use in a developing country like Bangladesh. 
The application of the techniques has been started since the last 
few decades, but due to some limitations, it is still backward and 
needs the immediate attention of the experts and government. In 
this paper we applied the RS and GIS techniques to study the 
land-use pattern change in the Bhanugach Reserved Forest and 
create a GIS data base for the study area. While working for the 
supervised classification, we faced problems with poor records in 
the respective beat office about their resources, which made our 
task troublesome. Considering this, the following recommenda¬ 
tions are made for the study area to which appropriate authority 


should pay necessary attention: 1) beat maps should be upgraded 
with proper features; 2) the maps have to be well organized so 
that those represent accurate features’ location and extent; 3) the 
database should be made up-to-date vegetation class-wise and 
their area context-wise as SPOT’s pixel size have a certain limi¬ 
tation of 20 m x 20 m; 4) plantations either pure or mixed should 
be defined with their age and density; 5) more keen monitoring 
of the whole area should be conducted. The most important mat¬ 
ter is to obtain further observations immediately, keeping in 
mind to get the field closer to the image scene. 
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